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1 Abstract 



We report the first isoelectronic differential force measurements between a Au-coated probe 
and two Au-coated films, made out of Au and Ge. These measurements, performed at submicron 
separations using soft microelectromechanical torsional oscillators, eliminate the need for a detailed 
understanding of the probe-film Casimir interaction. The observed differential signal is directly 
converted into limits on the parameters a and A which characterize Yukawa-like deviations from 
Newtonian gravity. We find a ~ 10 12 for A ~ 200 nm, an improvement of ~ 10 over previous limits. 

PACS numbers: 03.70.+k, 12.20.Ds, 12.20.Fv, 42.50.Lc 
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Although gravity was the first fundamental force to be understood, the quest to unify it 
with the other fundamental forces has remained elusive. One of the reasons is the apparent 
weakness of the gravitational interaction at small separations. Consequently, a significant 
number of experimental searches for new forces over ultra-short distances has been performed 
□□□□□□□□□ . They have been stimulated by at least three different — but related — 
motivations: 

(a) Some unification theories, incorporating n compact extra spatial dimensions with 
characteristic size R(n), predict deviations from Newtonian gravity over sub- mm scales 

n n n 

[10, llH, 112(| . Some extra-dimensional theories characterize the deviations by a Yukawa- 
modified potential V(r) 3 11, l3> 



V(r) = V N (r)[l + ae- r / x ] , (1) 

where Vjv(r) is the Newtonian gravitational potential for two point masses separated by a 
distance r 3> R(n), and a and A ~ R(n) are constants. Thus, the extra-dimensional theories 
provide a parameter R(n), which for n > 1 is naturally small, and an associated 
constant a n which is relatively poorly constrained. 

(b) String theory and other extensions to the Standard Model predict the existence of 
new light bosons such as dilatons, moduli, radions, cosmons, and bulk gauge bosons 
The exchange of these massive particles leads to corrections to gravity as in Eq. with 
a ^> 1 and A, related to the boson mass m by A = h/mc, as large as a few microns. Hence, 
the limits on a may provide useful guidance in narrowing down the myriad possible models 
linking physics at very high energy scales to the much lower energy Standard Model. 

(c) Theories in which these hypothetical new forces arise from an inverse-power potential, 
V oc r~ p . An extensive review of such inverse-power-law forces has been given recently in 
Ref. 

Existing limits on such forces, when parameterized as in Eq. are relatively weak for 
several reasons. Most significantly, if A is small then the effective interacting masses are 
themselves necessarily small, and background disturbances play a relatively more important 
role. Secondly, for experiments with typical separations ~ 1/im the dominant background 



arises from the Casimir force 







141 ]. which is not only relatively strong over the relevant 
distances, but is also somewhat difficult to completely characterize at the required level of 
precision jlj. Absent any alternative, limits on a = a(X) have been inferred by developing 



detailed theoretical models of the Casimir interaction yielding the limit a ~ 10 13 for 
A = 150 nm. Further improvements in experiments and in our understanding of the Casimir 
force will lead to stronger constraints on a (A), but it is better to sidestep theory altogether 
by carrying out what amounts to a "Casimir-less" measurement. 

In this Letter we report isoelectronic measurements (IET) 0, Q], where the Casimir 
background is subtracted at the outset, thus avoiding the necessity to model the Casimir 
force. IET exploit the essentially electronic nature of the Casimir force, whereas gravity and 
hypothetical forces couple to nucleons and electrons. Hence, vacuum fluctuations cannot 
account for any significant difference in the forces between a probe and two materials with 
identical electronic properties. 

A schematic of our set-up is shown in Fig. ^ We compare the force differences over 
two dissimilar materials, Au and Ge, which have been coated with a common layer of Au 
of thickness d p Au = 150nm > A p , where A p = 135 nm is the plasma wavelength for Au. The 
fractional difference of the Casimir force between two infinitely thick metallic plates and 
two plates of thickness d p Au is ~ e~ 47rd ^/ Ap ~ 10~ 6 1(|. In our experiment, it translates to 
a difference AFc ~ 10~ 17 N in the Casimir force between the Au coated sphere and the two 
sides of the Au/Ge composite sample. Hence any differential signal that the probe detects 
as it oscillates over the underlying Au and Ge (which provide a large mass density difference, 
PAu — PGe = 13.96 x 10 3 kg/m 3 ) must be due to an interaction via either gravity or some 
new hypothetical force. Thus by directly comparing the forces on the Au and Ge substrates 
a limit on a (A) can be obtained, without having to resort to a theory of the Casimir force 
for real materials. The expression for this hypothetical force difference is 



AF hyp {z) 
K s 

K rt 



-An 2 Ga\ 3 e- z/x RK s K p , 

[PAu ~ (PAu ~ PCr) e' d -- /A - (per ~ P s) e^-***^] 

(pAu-pce)e-^ +d ^/ x (l-e-**A)- 



(2) 



where G is the gravitational constant, R ~ 50/iin is the radius of the sphere, K s (K p ) is 
a term associated only with the layered structure of the sphere (plate), p s , pc r , PAu, and 
PGe are the densities of the sapphire sphere, a Cr layer (used to increase Au adhesion to 
sapphire), the Au layers and the Ge layer, respectively. Thicknesses di for the different 
materials are given in Fig. 



Since the hypothetical forces under study are weak a high sensitivity force measure- 
ment is required. A microelectromechanical torsional oscillator (MTO) with a soft spring 
k (high force resolution) and high quality factor (Q) satisfies the required demands. The 
MTO has a low coupling with the environment Q,!^ « ~ 10 _9 Nm/rad, and Q ~ 10 4 . 

The films deposited on the sphere and the MTO were characterized by atomic force 
microscopy (AFM). A typical AFM line-cut at the interface between the Au and Ge layers 
is shown in Fig.^Jx The observed ridge ( a valley in some samples) arises from the imperfect 
alignment of the mask when depositing the Au and Ge. The analysis of the AFM images 
indicates the granular character of the samples, showing a maximum height difference of 
22 nm. The average lateral dimension of the grains was in the 100-150 nm range, although 
grains as large as 500 nm were observed. The uncertainty in the position of the zero height 
level (with respect to which the roughness is measured Jj) is ~ 0.2 nm. 

The experimental arrangement and calibrations performed are very similar to those pre- 
viously used to determine Casimir forces 

HQ. 

A voltage was applied to the sphere to 
eliminate the residual electrostatic force caused by the difference in work functions between 
the Au layer on the MTO and the sphere. The angular displacement of the MTO, deter- 
mined by measuring the difference in capacitance between the MTO and the underlying 
electrodes, AC = C r i g ht — Ci e f t , yielded the force acting on it. The sensitivity in the angular 
deviation is 56 ~ 10 -9 rad/vHz. 

The force sensitivity is improved if the measurements are performed at resonance, u = u a . 
In this case, the minimum detectable force is dominated by thermal fluctuations, 5F(u ) ~ 
SFthermai = ^- / b \J AkIz bT / (u Q Q) , where ks is Boltzmann's constant. Consequently, it is 
necessary to measure the effect of AF hyp at resonance. AF hyp can be described as the product 
of a function of the separation e~ z ^ x and a function of the a;— coordinate, the difference in 
the average mass density. Hence, we induced a vertical oscillation on the MTO such that 
the separation between the MTO and the sphere changed as z m = z mo + 6zcos(u z t), with 
z mo 3> Sz. We simultaneously moved the MTO along a direction parallel to its axis, such that 
the effective mass density under the sphere was p e // = P + +P where = (pA u ±PGe)/2, 
and H(t) is a square- wave function with characteristic angular period T x = 2iru~ l . At t = 
the sphere is positioned over the MTO on the Au/Au half. The Casimir interaction leads 
to a shift of the resonance frequency of the MTO from its natural oscillation frequency f Q 
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to f r [17[. By selecting uj z + uj x = ui r , AFh yp has a Fourier component at uj r given by 

J~hyp( z mo, ^r) = AF hyp (z mo ) X — X I\ {^~^J ' (^) 

where I\ is a Bessel function of the second kind. Th yv (u r ) is the only Fourier component 
with a significant signal-to-noise ratio, even though no parts in the system are moving at f r . 
We selected f x = uj x /2tc = u r /14Q7r = f r /70 ~ 10 Hz. Consequently, f z = u z /2tc = 69 f x . 
The phases of the different signals were chosen to simultaneously cross through every 
t xr = 70/ f r (t xz = 69/ f z ) for the signals at f x and f r (f z )- The synthesized signal at f r 
was used as reference, see Fig. Ek- The amplitudes were adjusted to provide a peak-to-peak 
lateral displacement of D ~ 150/mi and a vertical amplitude that ranged between 5z ~ 
10 nm at the closest separations to 5z ~ 50 nm at z mo = 500 nm. 

For each resonant period of oscillation of the MTO, T r = l/f r , ten equally spaced 
data points J r d(ti) were acquired, ti = iT r /l0, % = 1, ■ • ■ , 10. Simultaneously, V(ti) = 
J^diti) cos u r ti and Q(ti) = J-a{U) smui r ti were determined. Averaging of J-'d, V, and Q was 
achieved by adding the signals for all different T r in the corresponding i — 1, ■ ■ • , 10 inter- 
vals. Furthermore, the summation over ti of V(ti) (Q(tj)) yields the in-phase T p (quadrature 
J-'q) Fourier component at uj r . 

A strong reduction of the random noise was achieved by increasing the integration times 
r. r was changed between 0.1 and 2000 s in a 1, 2, 5, 10, • • • sequence. Fig. El shows the 
results obtained at a separation z mo = 300 nm. Figs. |2t>-d show the observed behavior of 
J- p and T q in one sample for three characteristic values of r (the number of repetitions N 
decreased as r increased). Fig. |2jl, shows the behavior at r = 2, 000 s for all seven samples 
investigated. Several features should be noted in Fig. EJ (i) for each sample and all values 
of r \Tt\ = \j + T 2 q is constant within the statistical error. T v (J-'q) is the average of T v 
(J-'q) over the different repetitions; (ii) |JF T | remains constant within a factor of 2 for different 
samples; (Hi) the phase 6 = arctan(jF 9 /jF p ), however, assumes values close to either or 7r, 
indicating that the force is larger over either the Au or the Ge side of the composite sample, 
respectively. We assume that the forces measured correspond to either 9 = or 7r jis| . 

We calculated the force J-(r, z mo ) = J-t{t, z mo ) ± stv*(t, z^tp^* required to exclude, at 
the (3 = 95 % confidence level, any hypothetical force of the form given by Eq. (J2J). Here 
Sjv* (t, z mo ) is the mean square error of Tt(t-, z mo ), tp,N* is the Student's coefficient, and 
N* = 6 is the number of repetitions at r = 2000 s. The + (-) sign is is used for the in 



phase (-7T out of phase) cases. Fig. El shows the asymptotic behavior jF(r = 2, 000 s, z mo ) 
for all samples. \J-(r = 2, 000s, z mo )\ decreases as z mo increases, but not exponentially, as 
expected from Eq. (J3J). Independently of the origin of T{r = 2, 000 s, z mo ), however, the 
shaded region in Fig. El represents values of hypothetical forces excluded by our experiment. 

Although T{j = 2,000 s, z mo ) ^ 0, we do not believe it originates from new physics. 
As noted above, its z mo dependence is not exponential. More importantly, it changes sign 
depending on the sample under study, showing no correlation with the underlying Au and 
Ge layers. We therefore attribute its presence to the manifestation in the Casimir force 
of a change 8z' ~ 0.1 nm in z mo , i.e. z mo takes different values over the Au and Ge sides. 
With the AFM characterization, the height difference between the two sides is not known 
to better than 5z' ~ 0.1 nm [19j . This translates into a residual Casimir force difference 
| AFc(200nm)| ~ 15 fN, sufficient to explain the observed \J-(r, z mo )\. 

We can, however, rule out other sources for the observed background: (i) A motion of 
the sphere in a direction not parallel to the MTO's axis was ruled out by performing the 
experiments without crossing the interface, and at different z mo . The uncertainty 5y ~ 
3nm in the motion of the sphere over the D = 150 nm excursion yields a ~ 0.2 fN error. 
(ii) Local differences in roughness yield an estimate for the residual force to be £2fN, 

n 

while patch potentials [21] give a contribution £5fN. The effect of roughness and patch 
potentials can be further reduced by moving the sphere back and forth across the interface 
to points randomly selected, with dispersion large enough to average the local changes. This 
approach yielded the same results. (Hi) The effect of the ridge at the interface is also 
small. A Fourier analysis of its contribution yields $(200 nm) $3fN. As z mo increases, this 
contribution decreases rapidly, having a value $(500 nm) £0.1 fN. We note that for samples 
with a valley instead of a ridge the contribution is negligible, (iv) We also checked that the 
finite size of the sample does not affect the Casimir background or the Yukawa corrections. 
(v) Magnetic or gravitational (Newtonian) forces do not give a measurable background 
at uj t . The component of the magnetic force with a signal at u r is associated with the 
preferential presence of magnetic impurities in one of the sides of the engineered sample 
(either Ge or Au). This force, comparable to the magnetic interaction between isolated 
atoms, is much smaller than the sensitivity of our apparatus. The Newtonian gravitational 
attraction difference between the sphere and the composite Au/Ge sample is ~ 3 x 10~ 21 N, 
independent of separation. Hence, this force not only is too small to be detected, but it also 
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does not provide a Fourier component at LO r . 

The results shown in Fig. 3 can be used to obtain more stringent limits on hypothetical 
forces since these must be less than, or at most equal to, the observed background. The values 
of forces in the shaded region in Fig. |3]have been experimentally excluded. Furthermore, 
the 6 = and = n cases yield very similar absolute values for J 7 . Hence, for a given 
z mo , we use the smallest \T{r = 2,000s, z mo )\ as the maximum allowed hypothetical force. 
Using this value in the left-hand side of Eq. Q we obtain an a (A) curve. Repeating this 
procedure for different z mo we obtain a family of curves whose envelope provides the strictest 
limits arising from our experiments. This curve, together with previous results, is shown 
in Fig. which also shows the regions in A, a phase space where different models predict 
the existence of new forces |22j. Our realization of a "Casimir-less" IET yields a ~ 10-fold 
improvement (in the [30, 400] nm range) on existing limits for Yukawa-like corrections to 
Newtonian gravity. This has significant consequences for models of moduli exchange, as 
proposed by supersymmetry, by further constraining the supersymetric parameters. We 
believe that our direct, improved experimental test at submicron separations will continue 
to motivate theoretical development for this range of separations. We also note that our 
experiment can be improved by gluing the sphere to the MTO and oscillating the test masses 
over it. By judiciously designing the test masses, most of the problems associated with the 
background can be removed. In this scenario, limits on a down to 10 6 can be achieved at 
separations z mo ~ 100 nm. 
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across the interface and for stimulating discussions. R.S.D. acknowledges financial support 
from the Petroleum Research Foundation through ACS-PRF No. 37542-G. The work of 
E.F. is supported in part by the U.S. Department of Energy under Contract No. DE- 
AC02-76ER071428. 
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FIG. 1: (a) Scanning electron microscope image of the MTO with the composite sample deposited 
on it. Inset: Schematic of the sample deposited on the MTO. The coordinate system used in the 
paper is indicated, (b) Diagram of the experimental set-up. Inset: AFM profile of the sample across 
the interface. dc r = 1 nm ; d s Au = 200 nm: Au layer thickness on the sapphire probe; dc e = 200 nm: 
thickness of both the Ge (indicated in orange) and bottom Au layers on the substrate; dpt = 1 nm: 
thickness of the Pt film (used to avoid diffusion of Ge in Au) . There is also a dxi = 1 nm layer of 
Ti deposited to increase adhesion to the MTO. 

FIG. 2: (a) Integration time, r, dependence of the signal. Data at the shortest integration time 
have been scaled by A = 20. Data were obtained at z mo = 300 nm. The dotted line is the measured 
reference signal at uj r , adjusted to fit on the same scale, (b) through (d) Measured T v and T q 
signals with r = 2s (N = 300), 200s (N = 30), and 2,000 s (N = 6), respectively. Each point is a 
repetition of the experiment. The circles in (b) and (c) are centered at (J- p ,J- q ) and have a radius 
given by the calculated thermal noise over the relevant r. The measured noise, obtained at large 
z mo , exceeds the thermal noise by ~ 20%. (d) Results at r = 2000 s for all samples. 

FIG. 4: Values in the A, a phase space excluded by the experiment. The red curve represents the 
limit obtained in the current experiment. Curves 1 to 5 were obtained by Mohideen's group [3], 
our group Lamoreaux Kapitulnik's group [jj, and Price's group respectively. 



FIG. 3: Dependence of the peak-to-peak Fourier component Tt(j = 2,000s,2; mo ) on z mo for 
all measured samples. The shaded region in the figure represents values of hypothetical forces 
excluded by our experiments. Colors and symbols for the different samples are the same as those 
used in Fig. |2Jl. The upper (lower) hatched area is defined by the points with minimum absolute 
value of T(t, z mo ) = F T {r, z mo ) + s n *(t, z mo )tp, N * (F(t, z mo ) = T T {j, z mo ) - s n *(t, z mo )tp >N ,). 
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